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We investigate effects of finite core-spin in two-neutron halos in the three- 

o : 

body model consisting of two neutrons and a core. The states in the neutron- 
core subsystem can then have an additional splitting due to the coupling of 

•rH , 

the neutron and core spins. The connection between the structure of the to- 
tal system and the states of the neutron-core subsystem depends strongly on 
this spin splitting. However, the spatial structure is essentially unchanged. 
The neutron momentum distribution in a nuclear dissociation reaction will 
consist of a broad and a narrow distribution arising from the two spin split 
neutron-core virtual states. The large distance behavior is investigated and 
the conditions for the Efimov effect turn out to be more restrictive. The 
magnetic moment is equal to the core- value independent of spin splitting and 
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energies of the neutron-core states. The excited states of opposite parity 
depend strongly on the spin splitting. The parameters in our numerical ex- 
amples are chosen to be appropriate for n Li and 19 B. 
PACS numbers 21.45. +v, 21.60.Gx 



I. INTRODUCTION 



The existence of halo nuclei is by now established |1]-[3| as weakly bound and spatially 
extended systems. A rather accurate description is obtained by use of two- and three-body 
models. Criteria for the occurrence of halos and the related asymptotic large distance behav- 
ior were recently discussed [§-0]. These investigations assumed many simplifying approxi- 
mations, which allowed extraction of general properties and therefore provided insight into 
the general nature of such systems. Of particular interest are bound three-body systems 
where furthermore all two-particle subsystems are unbound. These so-called Borromean 
systems [Q appear for a substantial range of two-body potentials M. 

The nucleus 11 Li has played an especially dominating role in the developments. The ap- 
proximate description as a three-body system consisting of two neutrons outside the 9 Li-core 
was only recently suggested fpTol — 1 2 1 . Since then a series of improved calculations appeared, 



see for example || . The most detailed and accurate of these three-body computations are 



based on the Faddeev equations |fL3|| . In spite of the large number of publications the pos- 
sibly significant consequences of the finite core-spin s c has so far not been considered at all 
in three-body calculations. The neutron- and the core spins can couple to a total angular 
momentum of s c ± |, while the neutron-core motion still remains in the same relative orbital 
state. If the interaction is independent of these couplings, the structure corresponds to that 
of a spin-zero core. This is, however, highly improbable considering the well known appre- 
ciable spin dependence both of the strong interaction and of the observed bound states and 
resonances throughout the periodic table. The states in neighboring nuclei like 12 B, 12 N and 
8 Li as well as computations are in particular suggestive for a similar spin splitting in 10 Li, 



2 



see ref. [0,0 • 



The purpose of this paper is to investigate the consequences of finite core-spin for halo 
nuclei. The resulting spin splitting is probably an important ingredient for understanding 
the relation between the ground state structures of the neutron-core subsystem and the total 
three-body system, e.g 10 Li and n Li. Another effect is that the peculiar Efimov structure 
16| , which might show up in halo nuclei, is hindered by the spin splitting of the virtual s- 



states. As in our earlier papers |1|||7|], we want to establish the general features, qualitatively 
or quantitatively as best we can, and apply the results on realistic examples. We shall again 
for convenience assume that all three particles are inert and all the core degrees of freedom 
are therefore frozen. We shall furthermore assume throughout the paper that the total 
angular momentum J for the ground state is equal to s c and that the total parity equals 
that of the core. 

The paper is the third (after ||[7|]) in a series of papers discussing various aspects of three- 
body halos. After the introduction we sketch in section 2 the general theoretical method 
used to compute the properties of the three-body system. The interactions between the 
different two-body subsystems are parametrized in section 3. In section 4 the asymptotic 
form of the Faddeev equations is derived for interactions causing the spin splitting and 
in section 5 is discussed the large distance behavior in connection with the Efimov effect. 
The resulting spatial structure and the measured momentum distributions are discussed in 
section 6. Magnetic moment is investigated in section 7 and excited states of 1~ character 
on top of the ground state in section 8. Finally we give a summary and the conclusions in 
section 9. 

II. METHOD 

The Hamiltonian of the system, where the center of mass kinetic energy is subtracted, 
is given by 
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3 2 p2 3 

1=1 ' i>j=l 

where mj, and p« are mass, coordinate and momentum of the i'th particle, are the 
two-body potentials, P and M are respectively the total momentum and the total mass of 
the system. We shall use the Jacobi coordinates basically defined as the relative coordinates 
between two of the particles (x) and between their center of mass and the third particle 
(y). The precise definitions and the corresponding three sets of hyperspherical coordinates 
(p,a,Q x ,Q y ) are elsewhere defined, see e.g. 0-§f. Here p (= y/x' 2 + y 2 ) is the generalized 
radial coordinate and a, in the interval [0, tt/2], defines the relative size of x and y , Q x and 
Q y are the angles describing the directions of x and y. One of these sets of hyperspherical 
coordinates is in principle sufficient for a complete description. The volume element is given 
by p 5 dfldp where dfl = sin 2 a cos 2 adadQ x dQ y . 

The total wavefunction \1/ of the three-body system is written as a sum of three compo- 
nents each expressed in terms of one of the three different sets of Jacobi coordinates: 

* = X> (0 (x«y<). (2) 

i=l 

This three-component wavefunction is flexible and allows a description of different three- 
body structures by means of rather few angular momenta in each component. These wave- 



functions satisfy the three Faddeev equations [13 



(T - E)ipW + Vjk y,® + ipti) + ^W) = o, (3) 

where E is the total energy, T is the kinetic energy operator and {i, j, k} is a cyclic permu- 
tation of {1,2,3}. Any solution to eq. (||) is via eq. (0) also a solution to the Schrodinger 
equation, but the Faddeev equations, which in practice almost inevitably appear as integro- 
differential equations, are much better suited for descriptions of delicate correlations. 

Each component ip® is now for each p expanded in a complete set of generalized angular 
functions $W(p, f^). 

^ (i) -iE/«(*(^i). ( 4 ) 

P n 



where p _5//2 is the phase space factor. The angular functions must be chosen to describe 
correctly the asymptotic behavior of the wavefunction and in particular the behavior corre- 
sponding to possible bound two-body subchannels. This is essential for low energies where 
the correlations must be described very accurately. These functions are now chosen for each 
p as the eigenfunctions of the angular part of the Faddeev equations: 

~A 2 $« + V Jk (*® + *J> + $L fc >) = |^(p)^), (5) 

where {i,j, k} again is a cyclic permutation of {1, 2, 3} and A 2 is the p-independent part of 
the kinetic energy operator defined by 

m nn h 2 1 12 nn ^ ( <S/2 ^ 1 15\ 

Explicitly the generalized angular momentum operator A 2 is given by 

i o 2 i 2 i 2 

A 2 = — : — — - sin q. cos ex H — ; — % h — y - 4 (7) 

sm a cos a acr sin a cos^ a 

in terms of the angular momentum operators I 2 and l 2 related to the x and y degrees of 

freedom. This procedure is particularly convenient, since all angular variables are confined to 

finite parameter intervals and therefore correspond to discrete eigenvalue spectra. The idea, 

here applied to the Faddeev equations, was exploited previously to solve the Schrodinger 

equation in atomic physics [|l7j] and recently also in studies of the triton . 



The radial expansion coefficients f n (p) are component independent, since p is indepen- 
dent of which Jacobi coordinates are used. Insertion of eq. (f|) into eq. (|]) and use of eqs. (|J) 
(|) and (|7|) then lead to the following coupled set of "radial" differential equations: 
/ d 2 2mE 



V dp 2 h 2 

where the functions P and Q are the following angular integrals: 

d 



(8) 



Pnn'(p) = E jdm^(p,n)-^)( P ,n), (9) 



Qnn'(p) = E / ^^^fe^O^)- (10) 



<9p : 



Often very few terms of the expansion in eq. (f|) is needed to get sufficient accuracy in 
calculations of the lowest lying states. On the other hand the method is not suited for highly 
excited states where many high angular momenta are needed in the expansion. Systems 
with many crossings of the A-values as function of p are also numerically difficult to solve 
accurately by this method. 



III. TWO-BODY POTENTIALS 

Loosely bound systems are mainly sensitive to the low-energy properties of the potentials 
and we shall therefore use relatively simple potentials reproducing the available low-energy 
scattering data and still allowing extensive three-body calculations. The neutron-neutron 



interaction in the singlet s-wave is assumed to be the same as in ref. ||10|| , where the radial 
shape is a gaussian. In the extension to triplet p-waves we include spin-spin, spin-orbit and 
tensor terms. The resulting potential is then 

s nl " s n2 + VtS 12 + V so \ nn • S nn )exp(-(r/6 nn ) 2 ), (11) 

where s nl and s n2 are the spins of the neutrons, s nn = s nl + s n2 and Si 2 is the usual tensor 
operator. The four strength parameters and the range of the gaussian are then adjusted to 
reproduce the following four scattering lengths a and the s-wave effective range r e 



aCS ) = 18.8 fm, r^So) = 2.76 fm, 
a( 3 P ) = 3.6 fm, a( 3 Pi) = -2.0 fm, a( 3 P 2 ) = 0.30 fm. (12) 

The resulting values are given by 

Vc = 2.92 MeV, V ss = 45.22 MeV, V T = 26.85 MeV, V so = -12.08 MeV, 

b nn = 1.8 fm, (13) 



which by use in eq. (|12|) leads to the proper low-energy behavior both in the s-wave and the 
three different p-waves. 
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The neutron-core potential, which also is assumed to have a gaussian radial shape, is 
directly related to the resonances or virtual states in the subsystem. Their positions in 
energy are crucial in the computations and it is necessary to be able to adjust s-, p- and 
d-waves independently. We therefore use the parametrization 



■nc) )> 



= V s (l + 7s s c • s n )exp(-(r/b n 
V® = (Vt + V^s n • Uexp(-(r/&,„ c ) 2 ), (14) 

where / = 1,2 corresponds to either p- or d-waves, V s ,Vi,V$ , 7 S and 6 nc (=2.55 fm) are 
constants and s c , s n and l nc are the spin of the core, spin of the neutron and the relative 
neutron-core orbital angular momentum. The spin splitting is for simplicity only introduced 
in the s-state, but could as well be important for the other partial waves. Both the spin-orbit 
partners of higher orbital angular momenta may be independently adjusted. These effective 
potentials must take into account that the orbits occupied by the neutrons in the core are 
excluded by the Pauli principle. This is either done by choosing shallow potentials unable 
to hold bound states or by considering higher lying three-body states. 

We first choose V s = V p = -7.8 MeV, V® = and 7, = 0, which results in a total 



binding energy of 0.3 MeV corresponding to n Li. Maintaining this binding energy of the 
three-body system and varying V s and V so , we obtain the curve in fig. 1 relating the virtual 
si/2-state and the pi/2-resonance in the neutron-core system. When one of these states is 
very low the corresponding component is dominating in the three-body wavefunction. The 
straight lines between the calculated points are inserted to guide the eye. 

The spin splitting is now investigated by varying 7s and V s while keeping the binding 
energy at 0.3 MeV and V so at the values corresponding to the square and the vertical line in 
fig. 1. The resulting connection between the virtual s-states are shown in fig. 2 for these two 
cases. When the energy of one of the states increases the other decreases due to the inherent 
coupling in the Faddeev equations. A finite p-state admixture (the squares) increases the 
energies of the virtual s-states to maintain the total binding energy. In both cases the energy 
of one of the spin-split virtual s-states in the neutron-core system may occur anywhere below 



0.1 MeV. 

The angular momentum dependence is now investigated by using a <i 5 / 2 -resonance instead 
of the pi/2-resonance in the neutron-core system. This might be the relevant structure for 
19 B, where the 17 B-core probably has an essentially fully occupied p-shell. The total angular 
momentum of 17 B is 3/2 which also is the expected value for 19 B. This is consistent with 
l c = Inn = s nn = 0, which normally would be the state of lowest energy in our three-body 
calculation. The two-neutron separation energy of 19 B is estimated to be 500 ± 450 keV 



, see ref. ||20|| , and could in fact be anything below 1 MeV. Therefore we vary the energy 
along lines relating the neutron-core resonance energies by E(g?5/2)=xE(si/2), where x then 
is a number between and oo. Varying V s and Vd we obtain the total three-body energy 
shown in fig. 3 as function of E(s!/ 2 ) for several values of x. The p-waves are not included 
in the allowed configuration space. This is equivalent to a choice of p-strength V p , which 
places the p-states at high positive energies. The spin orbit strength Vffi = —10 MeV and 
the spin splitting parameter 7 S = are kept constant in the calculation. The starting point 
is V s = —7.87 MeV or V s = —7.40 MeV corresponding to a virtual s-state at 50 keV or 
100 keV, respectively. The total three-body energy decreases almost linearly with increasing 
resonance energies for all x- values and remains in all the cases considered within the limits 
of the estimate. 

The total binding energies for the black triangles in fig. 3 are around 0.45 MeV and 
0.75 MeV. The effect of spin splitting in the s-state is now investigated by increasing 7 S 
from zero while keeping the other parameters at values corresponding to these points. The 
strengths of the two potentials for a total spin of s c ± 1/2 are proportional to the related 
expectation values (1 + 7 s s c • s n ), which respectively is 1 + 7 s s c /2 or 1 — 7 s (s c + l)/2, see 



eq. ([14]) . The effective radial three-body potential is for small p a linear combination of these 
two two-body potentials with their respective statistical weights and it converges for large 
p towards the hyperspherical spectrum, i.e. values independent of the potential. This radial 
potential is consequently almost independent of j s and the resulting three-body energy is 
then also almost independent of 7 S . The connection between the energies of the two virtual 
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s-states corresponding to the spins s nc = 1, 2 of the neutron-core system is shown in fig. 4 as 
function of 7 S for interaction parameters corresponding to the black triangles in fig. 3. The 
qualitative behavior is similar to that of fig. 2, since the total energy remains approximately 
unchanged. 



As seen from eq. @ the angular eigenfunction A(p) is essentially an effective potential 
determined by eq. (|5|). The generalized angular functions $^(p, fij) is a sum of components 
with different angular momenta of the subsystems coupled to the total angular momentum. 
The radial equations at large p do not couple neither these components nor those belonging 
to different Jacobi coordinates with the important exception of the states with zero two- 
body orbital angular momentum. These s-states are often essential components in low lying 
states of interest and we furthermore only introduced the spin splitting in s-states. We shall 
therefore first concentrate on s-states where the angular dependence in eq. ([D is reduced to 
contain only a. 

The spin dependence of the i'th component of the total wave function in eq. @ is denoted 
xf\ where s is the intermediate spin obtained by coupling of the spins (sj and S&) of particles 
j and k. The spin s is afterwards coupled to Sj to give the spin of the core s c which in our 
case equals the total spin of the system. For each channel the angular part can then be 
written as 



IV. ANGULAR EIGENVALUE EQUATIONS FOR S-STATES 




(15) 



<5> {2 \p,a 2 ) 






(16) 



sin 2« 2 



sin 2«2 



$ (3) (P,«s) 




(p, as) 



X, 



(3) 

s c -l/2 



+ 




(P,a 3 ) 




(17) 



sin 2a 3 



sin 2a 3 



9 



where the channel i=l is defined by x between the two neutrons and where (j)^ is the spatial 
part of the wavefunction related to x^ ■ Since the total wave function must be antisymmetric 



(3) 

under interchanging of the two neutrons we have Ar ' 



Sc-l/2 



*> (2) and r/> (3) - -cb {2) 

? Sc -l/2 anCl V Sc +l/2 — ( Ps c +l/2- 



Substituting eqs . (|15|) — ([T7D into eq. (|5|) with i — 1,2 and subsequent multiplication from 

4 1} (p,ai), <j>s}-y*(p, "2) and <f>fj +1/2 {p, a 2 ): 



the left by Xo \ Xs}-x/2 anc ^ ^i!+i/2' we OD t arn the following set of equations for the functions 



d 2 



sin 2«i da 2 



• o w ^ 0o 1} (p, ai) 2 / ■ \ 
sm2ai - A(p) — — hp v nn (p sm a\) 



sin 2a 



sin 2ai 



+ 



(18) 



„ C-l/2^' a 2) ^+1/ 2 (P 5 «2) ^Ll/ 2 (P.«3) 
W : — „ h O2 : — ~ r Oi- 



t(2) 



sin 2a 2 



1 2 
sin 2a2 da 2 . 



sin 2a2 



sin 2a3 



+ C 2 



0+1/2 



(p,«3) 



sin 203 



0, 



^l/ 2 (^«2) 



^(p.ai) , C-l/2(^«2) ^_i /2 (P»«3) 



sin 2qg 



A (p) ) ' gc J!ol ~ +P 2 <x2 ) _ 1 /2l«nc(psina 2 )|x2 ) _ 1 /2> 



sin 2a 2 

(2) 



sin 2a 



+ 



sin 2a?2 



+ C 3 



sin 203 



C 



sin 2«3 



e> 2 



sin 2a2 Sa 2 



sin 2a2 — A(p) 

.(2) 



^+1/2 (P» "2) 

sin 2a 2 

(.(2) 



+ p 2 (Xsc+i 12 1 u «c(p sin a 2 )| X^+l/2 ) 

(2) 



0, 



,(2) 



(19) 



_ tfWl) ^ ^ c +l/2(P> a 2) ^ e -i /2 (P."3) ^ C+ 1/ 2 (P'«3) 

^2 — ; — 1 : — 04 ; — ^ ^3 



0, 



(20) 



sin2ai sin2a 2 sin2a3 sin2a3 

where Vij(x) = Vij(x/aij)2m/h 2 , ajk = [mjmk/(m(rrij + m^))] 1 ^ 2 , A(p) = A(p)+4 and where 
we used that the neutron- neutron interaction {v nn = V23) is spin independent for s- waves and 
that the neutron-core interaction (v nc = v\i = Vi 3 ) is diagonal in spin space. The coefficients 
Ci, i=l,2,3,4 depend on the spin of the core, and take the values: 







(2) 


c 2 = 


(xP 


(2) 


c 3 = 


(xt 


1/2IX 


c 4 = 


(xt 


i/ 2 lx 



(3) 



,(2) 



X s ) 



(21) 



Eqs.([l9|) and (p0| ) are identical to the ones obtained by substituting eqs.(|l5D-(|T7D into eq.( 

(3) (3) 

with i=3 and the subsequent multiplication from the left by X Sc -i/2 an< ^ X Sc +i/2- 
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The three spatial wavefunctions </>g and 4>^±i/ 2 are now f° r each p determined by 



eqs. (p!q)— (|20|) and the following boundary conditions 

tf \p, a = 0) = « = tt/2) = 02 } ±1/2 (p, a = 0) = <j> s % /2 (p, a = vr/2) = 0. (22) 

For p = the eqs.(|I"8|)-(|2"0|) clearly decouple and the solutions are 

<f>g]±i / 2 {p — 0, of) = 0q (p = 0, of) oc sin (2na), n=l,2, 3, ••• (23) 

and the eigenvalues are given by 

A = \(p = 0) = A(p = 0) -4 = AT(AT + 4), K = 2n - 2, (24) 

which is the well known hyperspherical spectrum. 

The behavior of A for small p-values can now be obtained in first order perturbation 
theory. The easiest is to add the three equations in eqs.(p^)-([20|) and thereby obtain the 



equivalent Schrodinger equation, which in turn then gives the following first order correction 
to A 



X(p) =X + p 2 

2s, + 2 



Vnn{p = 0) + - {xf}_ V2 Kc{p = 0)\xf^ 1/2 ) 



+ 



J (Xs c+ l/ 2 \ V nc(p = 0)\x K J +1/2 ) 



(25) 



2s c + 

where we used the spin independence of the neutron-neutron interaction. Eq. ( pS| ) reduces for 
vanishing core-spin s c = and spin independent interactions to the usual small p expression, 
where only the sum of the interactions enters, see ref . [0] . 

V. THE EFIMOV EFFECT 

The large distance behavior of the A- values reflect characteristic properties of the three- 
body system and its subsystems. In particular the necessary and sufficient condition for 
the occurrence of the Efimov states is simply that Aoo(= A(p — > oo)) is smaller than —4 or 
equivalently A < 0. Therefore eigenvalues and eigenfunctions for large p must be obtained 
from eqs.fll8D- (|20D and we shall follow the procedure in ref. Jl6| . 
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The wavefunctions depending on a k are first expressed in terms of the i'th set of hy- 
perspherical coordinates and the equation are subsequently integrated over the four angular 
variables Xfc and fi^. This amounts for s-waves to the substitution formally expressed by 
the operator R ik defined by 

1 1 rw/2-\ir/2-ipn,-ai\ 

Rik^=^-r= r • /o \ sm(2a k )ijj(p, a k )da k , (26) 

sm(2^ fc ) sm(2a i ) J\ Vih -ai\ 

where the transformation angle (fi k is given by the masses as 



Vik = arctao | (-l)-, p (m ' +m * +m3) , (27) 
1 V rriirrik J 

where (— l) p is the parity of the permutation {i,k,j}. The short ranges of the potentials 
confine a k to be smaller than the range tq of the potential divided by p. (The precise 
definition of r is given by p 2 v nc (r ) = A). The effect of the R-operator is then for large p 
approximately given as 

2a 

R ik^ ~ ^77^ ^(fi, <Pik) ■ ( 28 ) 

The wavefunction (p^ vanishes faster than 4^b}±i/2 P rov ided the scattering lengths a± 
related to the two relative spin states of the neutron-core system both are larger than that 
of the neutron-neutron system. We can therefore omit eq. (|18|) and neglect 0^ in eqs. ( |I"9D 
and (pOD and thereby for p ^> ro to first order in cto = ?"o/p obtain the simplified eigenvalue 
equations 

( d 2 



da 



Kp) ^-1/2 (P>«2) + P 2 (x ( s 2 ]_ 1/2 \vnc(psma 2 )\x { ^-i/2) 



^2-1/2 (P> "2) + j^ (C3< ^-l/2(/W) " C 4<I>?]+1/2(P> <P)) 



, (29) 



^ 2 - *(p) ) ^!+l/2^' Q!2 ) + ^ 2 ^2+ 1 / 2 l«nc(p sin a 2 ) |x2+ 1/2 > 



9a 



, (30) 

where = arctan</A c (yl c + 2) is given in terms of the nucleon number A c of the core. 
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We now proceed by dividing the a 2 -space into two regions I and II. In region I (ck 2 > 
a o — r o/p), where the potential is negligibly small, the solutions to the previous equations 
are: 

</Si/ 2 (P> "2) oc sin {{a 2 - |)V^) , (31) 

where we have imposed that the solution has to vanish at the boundary a 2 = ir/2. 

In region II (a 2 < cxo), where A/p 2 is negligibly small, we have two inhomogeneous 
differential equations with the solutions: 

«£-?/2(P. « 2 ) = «-(pa2) " (C 3 <3 /2 (P, ¥>) " C 4 ^ 1/2 (P, V)) , (32) 

^!'+f/2(P. «2) = n+(pa 2 ) - (-^^^^(p, ¥») - C 3 ^J 1/2 (p, , (33) 

where the inhomogeneous solutions, proportional to a 2 , easily are found due to the linear 
dependence of a 2 . The homogeneous solutions w±(pa: 2 ) are the wave functions describing 
the state of zero energy in their respective potentials (x^±i/ 2 |^nc(p a 2)|x^±i/ 2 ^' O^ide the 
range of the potential (pa 2 > r ) these wavefunctions take the form w±(pa: 2 ) oc pa 2 + a±, 
where a± are the scattering lengths of the two neutron-core potentials corresponding to the 
two different spin couplings. 

The eigenvalue equation for A now arises by matching of the derivative of the logarithm 
of the solutions at the division point a 2 = « = tq/ p. Keeping only the lowest order in r /p, 
and assuming that both scattering lengths are much larger than r (a± ^> r ), we finally 
obtain 

Acos 2 £y/i) + sin 2 + h*- + V^sin (7^) (34) 

Z lL — U-\- Z Z CL— Q-j- 

+7T^—r(— ~ — )— V sin ( J VA) sin (fa - J)\/A) ^- sin 2 (fa - ?-)y/i) = . 

2s c +l V a_ a+ v sin2<£ V 2 V ; vv ^ 2 ;v ; sin 2 2<p 2 ;v ; 

The dependence on the core-spin s c appears only through the scattering lengths a± and the 
parameter C 3 = — l/(2s c + 1). 

When one or both of the scattering lengths are finite and negative, which means that at 
least one bound state exists in the neutron-core subsystem, the lowest A-eigenvalue diverges 
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as —p 2 /a 2 , where a is the scattering length for the channel of strongest binding. When 
both a± are positive or infinitely large, we can distinguish between three different cases. In 
the first case both a+ and a_ are infinitely large. Then by keeping only the leading order, 
eq.(|34D reduces to 

A cos 2 VA) - -4tt sin2 (fa " ? ) V^) = . (35) 
2 sin 2(p 2 

which for n Li, where cp = arctan • 9, has the lowest solution at A = —0.00545. This 
value gives rise to an effective attractive potential in the radial equation eq.(|) given by 
— 0.25545/p 2 , which results in an infinite number of bounds states. This number is indepen- 
dent of the core-spin s c , and the only dependence on the nucleus is contained in the angle 
<p, see eq.(p7|). 

In the opposite limit, where p is large compared to both scattering lengths, i.e. a± <C p, 
the leading order approximation of eq.fl34]) is simply 



sin 2 (-VA) = (36) 

and the hyperspherical spectrum is recovered, i.e. A — 4 = K{K + 4), K = 0, 2, 4... 

The last case is characterized by one small and one large scattering length compared to 
p, i.e. < a_ C p < a + , or < a + < /J < a_. Now eq.(|34]) reduces to leading order to the 
form 



sin(-\/X) \/Acos (-VA) ± — — — sin ((</? - ^) y/X) 



2 



, (37) 



2 ' 2s c + lsin2^ " r 2' 

where the ± signs appear for a T < p <C a±. This equation is independent of a± and it has 
never any negative solutions. The lowest positive solutions for n Li are A = 0.702627 and 
A = 1.34930, respectively when a_ < < a + and a + p <^ a_. Since A always is positive 
this means that the barrier in the radial equation ((A — l/4)/p 2 ) never can give rise to the 
Efimov effect. This reflects the fact that infinitely many bound states only occur when at 
least two of the subsystems have infinitely large scattering lengths. If this extreme effect 
should be present both a± must be infinitely large. 
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The behavior of the A function is summarized in fig. 5, where the lowest A for n Li is 
shown. The scattering lengths have been chosen to be 10 2 fm and 10 fm, and the dashed 
(solid) line corresponds to a_ < a + (a + < a_). When both a_ and a + are much larger than 
p, both curves converge to the same value of A = —0.00545. In the opposite limit where p is 
larger than both scattering lengths, the curves both approach the value four corresponding 
to the lowest level of the hyperspherical spectrum. In the intermediate region, where p is 
much larger than one of the scattering lengths and much smaller than the other, A has a 
flat region. The constant value depends in this case only on which of the scattering lengths 
is larger than the other, but not on the specific values of a_ and a + . These A values are 
shown as dashed lines in the figure. 

It is important to emphasize that the conditions for the appearance of the Efimov states 
become more restrictive for a finite spin of the core. In principle, the effect appears when at 
least two of the binary subsystems have extremely large scattering lengths. We have seen 
above that both the possible couplings of the neutron-core spin states for halo nuclei must 
correspond to extremely large scattering lengths. If only one of the scattering lengths is very 
large the A function is not low enough to permit many bound states, i.e. occurrence of the 
Efimov effect. However, these conclusions do not affect halo nuclei with zero core spin. It 
should still be remembered that the most promising place in nuclei to look for the Efimov 
effect would be in a system, where one neutron is added to a pronounced one-neutron halo 
nucleus with zero core spin |7j]. 

VI. SPATIAL STRUCTURE 

The angular eigenvalues A are the effective radial potentials. The spectrum is therefore 
decisive for the spatial structure of the system. In fig. 6 is shown an example for interaction 
parameters corresponding to the square in fig. 1, i.e. realistic values for 11 Li. The beginning 
of the hyperspherical spectrum (K(K+4), K = 0, 2, 4...) is seen at both p = and at p = <x>. 
The lowest level decreases from zero at p = 0, goes through a minimum and returns back 
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to zero at infinity. The return to zero at large distances means that none of the two-body 
interactions is able to hold bound states. The negative pocket at smaller distances signifies 
an attraction in the total three-body system. The attraction is in this case able to hold the 
system in a bound state although no binary subsystem is bound. 

The angular structure corresponding to the eigenvalue spectrum is conveniently dis- 
cussed in the Jacobi coordinate set where x is the vector between the two neutrons. The 
components of the wavefunction for p = and p = oo corresponds to (l nn , l c , L, s nn , s c , S) 
= (0,0,0,0,3/2,3/2) for the lowest level. The five fold degeneracy of the K = 2 levels 
corresponds to the angular components (l nn ,l c , L, s nn , s c , S) = (0,0,0,0,3/2,3/2) with one 
node, (1, 1, 0, 1, 3/2, 3/2) and the three levels (1, 1,1,1, 3/2, S) with S = 1/2, 3/2, 5/2. 

The usual measure of the size of the system is the mean square radius, which for A c 
nucleons in the core is given by 



where R 2 RMS and R 2 RMS {coie) are mean square radii of the respective systems and the mean 
square of the hyperradius is given by < p 2 >= A c < r 2 > +2 < r 2 >, where r c and r n are 
core- and neutron-coordinates measured from the center of mass of the total system. 

In fig. 7 we show the mean square radii along the curves in figs. 2 and 4 as functions of 
energy of one of the virtual s-states. Clearly this measure of the extension of the system is 
essentially independent of the spin splitting in all cases. This is due to the fact that the spin 
splitting gives rise to very small changes of the angular eigenvalue spectrum. The lowest 
A functions are indistinguishable at large distances and only small differences occur at the 
smaller distances, where the wavefunction has small contributions. 

However, the size depends both on the binding energy and the amount of higher angular 
momentum states in the neutron-core motion. For a given binding energy, the energy of 
the virtual s-state decreases with increasing s-state content in the neutron-core subsystem, 
and then the spatial extension of the system increases due to the smaller centrifugal barrier. 
The essential parameter is the amount of s-wave in the relative neutron-core motion. The 




(38) 
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variation of the root mean square radius of n Li and 19 B as functions of the s-wave probability 
along is shown in fig. 8. The white circles represent the variation of the radius for n Li along 
the curve in fig. 1. The black circles represent the same variation for 19 B, but along the 
horizontal dashed line in fig. 3, which corresponds to a constant energy of around 0.65 MeV. 
The resulting almost linear dependence for n Li is varying within limits consistent with both 
computations |§ and the recent interpretation (Rrms( u ^) = 3.1 ±0.3 fm) of the measured 



reaction cross sections 21 



In the ground state the ratio between the weights of the spin split s-waves are predefined 
by the symmetry constraint in such a way that the parametrization in eq. flHD makes the 
A values nearly independent of 7 S . The resulting spatial three-body structure is therefore 
independent of the spin splitting provided the total energy and thereby the total s-state 
occupation probability remains unchanged. Any observable for the three-body system de- 
pending solely on the spatial structure is therefore unable to determine the spin splitting 
of the subsystem. Direct measurements of the neutron-core virtual states would of course 
in principle be able to give the information. However, angular momentum determination of 
these unstable nuclear states are needed. 

A recent suggestion to determine the s-state admixture in the three-body ground state 



of n Li is highly relevant in this connection [p2 ]. The idea is that in a nuclear break- 
up reaction the final-state interaction between the core and the neutron is decisive for the 
observed neutron momentum distribution. Low lying virtual s-states or equivalently strongly 
attractive s-state potentials will via the final state interaction lead to a narrow momentum 
distribution [p3fl . A higher lying s-state or a p-state with repulsive centrifugal barrier can 
only weakly influence the outgoing neutron. The necessary analyses require first of all the 
probability of finding the remaining neutron-core system in the low lying s-state. For a 
non-negligible spin splitting this probability equals the probability of finding two occupied 
s-states multiplied by the probability of survival (roughly the statistical weights independent 
of the spin splitting) of the low lying state. Therefore information about the spin splitting 
may be obtained from these reactions. 
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VII. MAGNETIC DIPOLE MOMENT 



Structural changes due to the spin splitting are most likely to show up in spin dependent 
observables. The natural first choice to study is the magnetic dipole moment, which in our 
three-body model is given by 

A = (j M = J\(g n s nn + g c s c + Z c —^—\ c ) z \JM = J), (39) 

flN A c + I 

where /i^r is the nuclear magneton, Z c is the charge of the core, l c is the orbital angular 
momentum of the core relative to the center of mass of the neutrons, J and M are the total 
angular momentum and its projection. The gyromagnetic factors are taken as g n = —3.8263 
for the neutron, g c = 3.4391/s c for n Li and the Schmidt value g c = 3.7928/s c for 19 B. The 
non-vanishing matrix elements of the operator in eq. fl39|) are then 



{iJnrJ'c)L(yS nn S^)S ', J J\ (s nn ) 2 \(l nn l c ^L( y S nn S c ^S , J J) 



'-iy+L+s+s>+s nn+ s c Jj( 2J + 1)(2S + 1)(2S , + i) Snn ( Snn + i)( 2Snn + i)/(j + i) 

(40) 



<7 1 <?' s Is 



J L J J I S' s c S 

where L is the total orbital angular momentum, l nn is the orbital angular momentum of 
the two neutrons, S (S') is the total spin. The expression for the other two terms of the 
operator in eq. ( j39|) can be obtained respectively by exchange of s c and s nn , and by additional 
exchanges of all spins and orbital angular momenta on the right hand side of the equation. 

The results of these formulae are given in table 1 for both our examples of 11 Li and 
19 B. The major contributions consistent with parity and angular momentum coupling rules 
are listed in the neutron-neutron Jacobi set of coordinates. All diagonal p-wave matrix 
elements except one are smaller than the s-wave matrix element. Since all these components 
are coupled in the three-body system one could expect that p-wave admixture will reduce 
the magnetic moment compared to the core-value in contradiction with measurements for 
n Li |24j|. However, the diagonal elements alone are very misleading due to the strong off 



diagonal matrix elements given at the bottom of the table. 



The result of direct numerical calculations for both n Li and 19 B show that the core- 
value of the magnetic moment always is recovered and it is therefore insensitive both to 
the spin splitting and to the amount of p- or d-wave admixture. The reason is that the 
basic neutron-core components in the present three-body model for 11 Li are the (si/2) 2 and 
(P1/2) 2 two-particle states. Both are fully occupied j-shells of total angular momentum zero. 
Their contributions to the magnetic moment therefore must vanish and the off diagonal 
elements are essential in this connection. We are therefore left with the core-value. Our 
model also allows other configurations, which however turn out to be essentially unoccupied 
in the final wavefunction. For 19 B, where the p-states are substituted by d-states, we also 
recover the core- value for the magnetic moment. The reason is that the new configuration 
(d5/2) 2 almost exclusively couples to the total spin zero and therefore also in this case only 
contributes insignificantly. Other components are apparently energetically too expensive. 

The accuracy of the strict three-body model is relying on the assumption of an inert core. 
Thus core degrees of freedom must be included when improvements beyond the three-body 
model are requested. The most obvious modification for 11 Li is to allow an admixture of the 
only particle stable excited state, an angular momentum 1/2-state [II], into the 3/2- ground 



state of 9 Li. Since all matrix elements of the magnetic dipole operator in eq. ( |3"9"D vanish 
between states of different core-spin, only the contributions in table 2 are able to modify our 
previous results. Although we assumed an unchanged gyromagnetic factor for the core, it is 
rather obvious that any s c = 1/2 admixture would further reduce the magnetic moment. A 
value of s c larger than 3/2 is needed, but not available in the discrete spectrum of 9 Li. The 
effect of higher angular momentum neutron-core states and a j-shell only partially occupied 
by the halo neutrons does not alter the conclusions above as seen from our example of 19 B. 

VIII. EXCITED STATES OF ELECTRIC DIPOLE CHARACTER 

The spectrum of excited states in halo nuclei is almost by definition limited to at most 
a few discrete states, but the continuum structure including resonances might also be in- 
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teresting. The 1~ type of excitations, continuous or discrete, are for such nuclei expected 
to differ significantly from corresponding states in normal nuclei. The only excited state in 
n Be has the largest measured B (El)- value to the ground state [p5[| . In general the dipole 
strength function is peaked at low excitation energies and it plays a dominant role in the 
analysis of Coulomb dissociation cross sections |I|,^,P7[] . The peak has even been suggested 
to be a new type of low lying excitation mode for halo nuclei [P8| , p^| . In any case it can be 
expected that these 1~ excitations are the lowest lying excitations in halo nuclei. 

When the total angular momentum as well as the spin of the core both are zero the major 
components in the wavefunction of the excited 1~— state are most conveniently described in 
the neutron-neutron coordinate system. The angular momentum 1 and the negative parity 
can be made by (l nn , l c , L, s nn , s c ) = (0, 1, 1, 0, 0) and (1, 0, 1, 1, 0) if we for simplicity exclude 
d- and higher partial waves. The finite core spin now dictates the additional coupling of s c 
to a total angular momentum of J = s c , s c ± 1. The parity is still the opposite of the ground 
state. 

For a core spin of 3/2 corresponding to n Li the total angular momenta are J = 1/2, 3/2 
and 5/2. They are shown in fig. 9 as functions of the spin splitting parameter j g . The ground 
state energy is also shown for comparison and as expected seen to be almost independent of 
the spin splitting. In contrast the energies of the excited states vary with j 8 . An increase 
is seen for J = 1/2, 3/2 and a decrease for J = 5/2. Crudely speaking the l~-state consists 
of one neutron in an s-state and the other neutron in a p-state relative to the core. By 
selecting the lowest lying of the (spin split) s-states the energy of the three-body state could 
therefore be expected to decrease with the energy of one of the s-states in the neutron-core 
subsystem, i.e. as soon as 7 S ^ 0. However, the 1~ states have admixtures, determined by 
a complex interplay of several factors, of spin-split s-waves different from the ground state. 
Due to the antisymmetry requirement it is impossible to construct a wavefunction with given 
total J and one neutron in the lowest neutron-core spin split s-state and the other neutron 
in a p-state. Both virtual s-states are necessarily present. The J=5/2 state has an excess of 
neutron-core states with angular momentum s nc = s c + 1/2 = 2, the J=l/2 has analogously 
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an excess of s nc = s c — 1/2 = 1 and the J=3/2 state is somewhere in between. 

From the interaction in eq. ( |T4| ) we see that a negative value of V^g implies that the 
energy of the neutron-core state with angular momentum s c + 1/2 is lower than that of 
s c — 1/2. In our specific case of s c = 3/2 this corresponds for positive 7 S (V s < 0) to a 
lower energy for the state of 2 than for that of 1. Thus the energy of the state with the 
largest content of 2 must exhibit the fastest decrease with 7 S . Therefore with increasing 7 S 
the J=5/2 state necessarily becomes lowest in energy followed by the J=3/2 state and with 
J=l/2 as the highest. The details of the resulting energies are seen in fig. 9. An increase of 
the ground state energy would bring the excited spectrum into the continuum, but a similar 
structure would remain. 

IX. SUMMARY AND CONCLUSION 

We discuss in this paper the effects of a finite core spin on the structure of two-neutron 
halo nuclei. In the previous two papers in this series we established the general connections 
between size and binding energy for loosely bound three-body systems. The possible asymp- 
totic (spatially extended and small energy) structures were characterized by the interactions 
between the constituent binary subsystems. The largest of the "particles" called the core is 
assumed to be inert and spinless and the corresponding degrees of freedom are frozen. When 
the core has a finite spin the two-body interaction between the core and another particle 
(typically a neutron) depends on the total spin of the two-body system. More complicated 
spin structures therefore seems to be possible even when the core still remains inert. The 
present paper, number III in the series, is devoted to study effects of such spin couplings on 
the three-body system. 

The recently developed method to solve the Faddeev equations is well suited for our 
purpose. We therefore first described the necessary ingredience of our procedure. The 
definitions and notation are then established in general and we proceed to consider the 
crucial two-body interactions. The neutron-neutron potential is needed in both relative s- 
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and p-waves and since the spin structure is under investigation we also have to distinguish 
between the singlet and triplet states with the different total angular momenta. Low energy- 
properties are still expected to be decisive and we construct therefore an interaction which 
reproduces the four phenomenologically extracted s- and p-wave scattering lengths and in 
addition the s-wave effective range. This is achieved by use of spin-spin, spin-orbit and 
tensor terms with a common gaussian radial dependence. The resulting effective low energy 
neutron-neutron interaction might also be useful in other connections. 

The unknown neutron-core effective potential is also parametrized with a gaussian radial 
shape. The minimum spin dependence include a spin-orbit interaction and also the central 
interaction is allowed to be different in each of the relative orbital angular momentum 
channels. We furthermore for s-waves include a term to distinguish between the possible 
couplings of the core- and neutron spins and we choose the scalar product of these spins. The 
resulting interaction is now flexible enough to study the influence on the three-body system 
of the spin splitting and energy dependence of the resonances in the binary subsystems. 

Equipped with an efficient method and the input properly parametrized we then explic- 
itly write down the angular eigenvalue equations when only s-waves are included. These 
states are often the most essential components in the solutions. We find analytically both 
small and large distance behavior of the angular eigenvalues. The conditions for the Efimov 
effect are simultaneously obtained and investigated in terms of scattering lengths for the two 
different spin couplings of the neutron-core subsystem. The extreme Efimov effect, where 
infinitely many bound states are present in the three-body system only occurs when both 
scattering lengths are infinitely large. Already the first of these excited states is very ex- 
tended in space compared to the ranges of the potentials and the conditions for its presence 
is therefore interesting even when no other excited state exists. These conditions may conve- 
niently be formulated in terms of bound state and virtual s-state energies of the neutron-core 
subsystem. Both the two spin couplings of s c ± 1/2 must correspond to energies between 
about —2 keV and 5 meV. Since the typical spin splitting of such s-states is fractions of an 
MeV, the occurrence of the first Efimov state is extremely unlikely for a halo nucleus with 
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a finite core spin. 

The spatial structure of the three-body ground state turns out to be essentially inde- 
pendent of the spin splitting of the virtual s-state. The reason is that either none or both 
s-states are needed in one specific combination, if a total spin zero state has to be formed 
from the spins of the two neutrons. Any other combination has the wrong parity. On the 
other hand the spatial extension depends on the amount of s-state admixture as noted in 
previous publications and easily understood from the general behavior of states with higher 
angular momenta or higher hyperspherical quantum numbers. 

The energy of the three-body ground state depends sensitively on the energies of the 
(spin split) virtual neutron-core s-states, which constitute the components of the total sys- 
tem. The connection between the properties of the neutron-core system and the total three- 
body system depends strongly on the spin splitting. This may be exemplified by 10 Li and 
11 Li, where a given s-state content in the 11 Li-wavefunction requires a specific combination 
the two virtual s-states in 10 Li. If one is low lying the other energy must be high in or- 
der to compensate and lead to the observed total three-body binding energy. Thus, both 
virtual s-states in 10 Li must be known, if the 11 Li-binding energy should be accounted for 
unambiguously. 

The neutron-core s-state admixture in n Li is, as recently suggested, possible to extract 
from the measured momentum distributions in nuclear break-up reactions. In the sudden 
approximation, where one of the neutrons is removed instantaneously, there is a finite prob- 
ability of finding the other neutron in the lowest lying virtual s-state. Due to the strong 
attraction this part of the wavefunction leads to a narrow observable momentum distribu- 
tion. 

The most obvious spin dependent observable is the magnetic dipole moment, which for 
11 Li is found to be fairly close to the core value measured for 9 Li. The results of our strict 
three-body calculations are that the core-value is reproduced independent of spin splitting 
and s-state admixture. The reason is that the spins of the two neutrons always couple to 
zero even when other possibilities are allowed. 
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The discrete or continuous spectrum of excited states of l~character on top of the ground 
state is of special interest for halo nuclei. One of the main components in the three-body 
wavefunction consists of one neutron in an s-state and the other neutron in a p-state relative 
to the core. Coupling of these angular momenta to 1 unit of h corresponds to an electric 
dipole excitation. More than one electric dipole excitation corresponding to different angular 
momenta is possible for finite core spin. The related hyperfine splitting of the El-excitation 
leads to predictable structure in the continuum. Analyses of for example Coulomb dissoci- 
ation cross sections may elucidate the structure of this 1~ excitation spectrum. 

The excitation energy of such states may for a given total energy depend strongly on the 
spin splitting of the virtual s-states. The reason is that the excited states are composed of 
mixtures of spin split s-states different from that of the ground state. On average the lowest 
s-state is preferred and the statistically weighted excitation energies therefore decrease when 
the energy of the lowest virtual s-state approaches zero. The absolute value of this average 
depends on the energy of the p-state resonance. 

In conclusion, we have investigated effects of finite core spin in halo nuclei. We found 
essentially unchanged spatial structure and magnetic dipole moments. Efimov states seems 
to be extremely unlikely in halo nuclei with finite core spin. The electric dipole excitations 
may depend rather strongly on the relative positions of the resonances and virtual states in 
the neutron-core subsystem. The spin splitting of the s-states in the neutron-core system 
is decisive to understand the proper relation between the total three-body system and its 
subsystems. Our numerical examples concentrated on 11 Li (neutron-neutron- 9 Li) and 19 B 
(neutron-neutron- 17 B). These examples are used for general illustration, but they are useful 
in their own right as model computations for these nuclei. 
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TABLES 



Table 1 

Diagonal contributions to the magnetic moment from two-body s- and p-states in the 
neutron-neutron Jacobi coordinate system. The total and core angular momenta and 
parities are given by J 71 " = 3/2~ and si = 3/2~. The non- vanishing off diagonal matrix 

elements are shown at the bottom of the table. 
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Table 2 



Diagonal contributions to the magnetic moment from two-body s- and p-states in the 
neutron-neutron Jacobi coordinate system. The total and core angular momenta and 
parities are given by J 71 " = 3/2~ and si = l/2~. The non- vanishing off diagonal matrix 
elements are shown at the bottom of the table. The gyromagnetic factors of the core are 
respectively g c = 3.4391/1.5 and g c = 3.7928/1.5. 



no 


*"nn 


k 


L 




S 


((Snn)z) 


<(Sc)*> 


((l c ) 2 ) 


M( Li) 


/I Ot-v \ 

M 19 b) 


1 





2 


2 





1/2 





-3/10 


9/5 


0.29 


0.19 


o 
L 


Z 


u 


Z 


n 
u 


1 /9 
1/Z 


n 
u 


— O/ 1U 


n 
u 


— u.oy 


n 7f; 

— V. 10 


3 


1 


1 


1 


1 


1/2 


2/3 


-1/6 


1/2 


-2.66 


-2.71 


4 


1 


1 


2 


1 


1/2 


-2/5 


1/10 


9/10 


2.25 


2.26 


5 


1 


1 





1 


3/2 


1 


1/2 





-2.68 


-2.56 


6 


1 


1 


1 


1 


3/2 


11/15 


11/30 


1/5 


-1.86 


-1.77 


7 


1 


1 


2 


1 


3/2 


1/5 


1/10 


3/5 


-0.21 


-0.20 






<3| 




|4> 










1/^/20 


0.12 


0.12 






<3| 




|6) 




-2/V45 


-2/V45 





0.46 


0.39 






<4| 




17) 




-2/5 


-2/5 





0.61 


0.52 






<5| 


A* 


|6) 










v/275 


0.34 


0.33 






<6| 


A* 


17) 










2/5 


12/55 


4/19 



28 



Figure Captions 

Figure 1 The relation between the virtual Si/ 2 -state and the pi/ 2 -resonance for a total bind- 
ing energy of 0.30 MeV. The curve is obtained by varying V s and while keeping 
V p = —7.8 MeV and 7 S = 0. The square, the diamond and the dashed lines (cor- 
responding to the asymptotic values) indicate points of interest in the following 
figures. 

Figure 2 The relation between the two virtual s-states (total neutron-core spin S nc = 1, 2) 
for a total binding energy of 0.30 MeV for 11 Li. The curves are obtained by 
varying V s and 7 S while the other interaction parameters are defined by the points 
indicated in figure 1 by the vertical line (circles) and the square (squares). 

Figure 3 The total energy of 19 B as function of the energy E(s 1 / 2 ) of the second virtual 
s-state. The energy of the d 5 / 2 -resonance in the neutron-core system is given by 
E((i5/2)=xE(s 1 / 2 ) where x is constant along the curves. The parameters in the 
calculations are V$ = -10 MeV, and 7, = while V s and V d are varied. 

Figure 4 The same as fig. 2 for 19 B. The parameters correspond to the black triangles in 
fig. 3 with the total energy equal to 0.45 MeV and 0.75 MeV. (The curves are 
obtained by varying 7,,.) 

Figure 5 Large distance behavior of the angular eigenvalue A as function of p for 11 Li. 

The scattering lengths (marked by the vertical arrows) are a_ = 10 2 fm and 
a + = 10 7 fm for the dashed line and a + = 10 2 fm and a_ = 10 7 fm for the 
solid line. The horizontal dashed lines correspond to the various limiting values 
described in the text. 

Figure 6 The angular eigenvalues as function of the radial coordinate p for interaction 
parameters corresponding to the diamond on fig. 1, i.e n Li. Only antisymmetric 
states in the neutron-neutron exchange are exhibited. 
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Figure 7 The root mean square radius of n Li and 19 B as functions of the energy of one 
of the virtual s-states for interaction parameters corresponding to the curves in 
figs. 2 and 4. The squares, circles and triangles have the same meaning as in those 
figures. 

Figure 8 The root mean square radii of n Li (open circles) and 19 B (black circles) as func- 
tions of the s-state admixture in the neutron-core channel. The core radii are 
-Ri?Ms( 9 Li) = 2.32 fm and -Rrms( 17 B) = 2.65 fm respectively. The interaction 
parameters correspond to the curve in fig. 1 and the horizontal dashed line in 
fig 3. 

Figure 9 The ground state and the excited l~-levels as functions of the spin splitting 
parameter 7 S . The interaction parameters are V s = —8.5, V p = —7.8 and 
Vl p) = 34.3 MeV. 
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